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Abstract 

In this paper we calculate the effects produced by the temperature 
in the renormalized vacuum expectation value of the square of the 
massless scalar field in the pointlike global monopole spacetime. In 
order to develop this calculation, we had to construct the Euclidean 
thermal Green function associated with this field in this background. 
We also calculate the high-temperature limit for the thermal average 
the zero-zero component of the energy-momentum tensor. 
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1. Introduction 



It is well known that different types of topological defects may have been 
created in the early Universe after the Planck time by the vacuum phase 
transition Ijl], 0. These include domain walls, cosmic strings and monopoles. 
Among them cosmic string and monopole seem to be the best candidates to 
be detected. 

A global monopole is a heavy object formed in the phase transition of 
a system composed by a self-coupling scalar field triplet ip"', whose original 
global 0(3) symmetry is spontaneously broken to U(l). The simplest model 
which gives rise to a global monopole, was presented by Barriola and Vilenkin 
[|]. The gravitational effects produced by this object may be approximated 
by a solid angle deficit in the (3+l)-dimensional spacetime whose line element 
is given by 

ds^ = -a^dt^ + \dr'^ + r'^{e + sin^ Od^"^) . (1) 
a 

Here the parameter = 1 — SirGrjQ is smaller than unity and depends on 
the energy scale 7]q where the global symmetry is spontaneously broken. The 
energy-momentum tensor of this monopole has a diagonal form and reads: 
TO = = (a2 - l)/r2 and T| = T| = 0. 

The non-trivial topology of this spacetime implies that the renormalized 
vacuum expectation value (VEV) of the energy- momentum tensor, {T^^{x))ren 
associated with an arbitrary collection of conformal massless quantum fields 
should not vanish [Q. The explicit calculations for the massless scalar ||^ and 
fermionic[^ fields have been already obtained. 

In the framework of the quantum field theory at finite temperature, the 
fundamental quantity is the thermal Green function, (^/^(x, x'). For the scalar 
field it should be periodic in the imaginary time with period /5, which is 
proportional to the inverse of the temperature. Because we are interested to 
obtain the thermal Green function, it is convenient to work in the Euclidean 
analytic continuation of the Green function performing a Wick rotation. So, 
we shall work on the Euclidean version of the monopole metric, 

ds^ = dr^ + ^dr^ + r\e + sin^ Odif^) , (2) 
a'' 

where we have absorbed the a-parameter redefining the Euclidean temporal 
coordinate. 
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The main objectives of this paper are to study the effects produced by the 
temperature on the renormahzed VEV of the square of the massless scalar 
field operator and the respective energy-momentum tensor in the spacetime 
defined by (0). In order to do that, we calculate first the respective Euclidean 
thermal Green function in this manifold, which is the most relevant quantity 
for the obtaining of these expressions. 

This paper is organized as follows. In section 2, we obtain the Euclidean 
thermal Green function Gj3{x,x') adopting the imaginary-time approach]^], 
using the Schwinger-De Witt formalism for a massless scalar field. In section 
3, we calculate the renormahzed thermal average value (0^(a:))/3. Because this 
term cannot be written in a closed form, its dependence on the temperature 
is not evident. So, in order to obtain some quantitative information about its 
behavior, we have to proceed a numerical evaluation. In section 4, we present 
a formal expression for the thermal average of the zero-zero component of 
the energy-momentum tensor (Too(x))^. Again, because of its complicated 
dependence on the temperature, only numerical calculations enable us to 
present its appropriate behavior. In section 5, we present our conclusions 
and remarks about this paper. Finally, we left for the Appendix some details 
of our numerical calculations. 

2. The Euclidean Thermal Green Function 

The Euclidean Green function associated with a massless scalar field in 
the global monopole spacetime has been obtained by Mazzitelli and Lousto 
a few years ago[^]. There this Green function was obtained assuming a non- 
minimal coupling between the field and the geometry through the scalar 
curvature, R = 2(1 — a^)/2r^, of this spacetime. 

In this section we extend this result to obtain the Euclidean Green func- 
tion at finite temperature in the monopole spacetime defined by (H). The 
thermal Green function, Gi3{x,x'), must obey the non-homogeneous Klein- 
Gordon differential equation, and be periodic in the "Euclidean" time r with 
a period f3 = I/k^T, kb being the Boltzmann constant and T the abso- 
lute temperature. Considering the non-minimal coupling, the scalar Green 
function must obey the differential equation: 

{n-^R)Gp{x,x') = -6'^'\x,x'), (3) 
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where □ denotes the covariant d'Alambertian operator in the metric defined 
by (0), ^ is an arbitrary couphng constant and d^^\x,x') is the bidensity 
Dirac distribution. 

Because our spacetime is an ultrastatic oneQ the Euchdean thermal Green 
function can be obtained by the Schwinger-De Witt formahsm as follows: 



Gf3{x,x') = / dsKp{x,x'] s) , (4) 
where the thermal heat kernel, Kp{x^x'\ s), obeys the equation 

(^^-□ + ei?)i^^(x,x',s) = (s>0), (5) 

is subject to the boudary condition 

limi^^(a;,x';s) = 5'^^\x,x'), (6) 

and periodic in the Euclidean time with period j3. 

Following the prescription given in the papers by Braden^] and PageH], 
the thermal heat kernel can be expressed in terms of the sum 

oo 

Kp{x,x]s)= ^ Koo{x,x' -nX[5]s) , (7) 

n=— oo 

where A is the "Euclidean" time unit vector. The zero temperature heat 
kernel, K^{x, x'; s), in this spacetime, can be factorized as 

Koo{x,x';s) = K(^i){t,t']s)K(^3){x,x;s) , (8) 

where i^(i)(r, r'; s) and K(^3){x, x'; s) obey, respectively, the differential equa- 
tions: 

(§;-^)'<ir,(ry:s)-0 (9) 

and 

~ + ^) = • (10) 

^An ultrastatic spacetime admits a globally defined coordinate system in which the 
components of the metric tensor are time independent and the conditions goo — 1 ^nd 
goi = hold. 
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Here the covariant spatial Laplace operator has the form 

is the square of the flat-space angular momentum operator. 



Also, as was poited out in |T^, the heat kernel can be expressed in terms 



of the eingenfunctions of the operators — □ + ^i? as 

Kix, x'- s) = Y. 0.(x)0:(x')e-^'^' , (11) 

(7^ is the corresponding positively defined eigenvalues. Writing 

+ = a20,(x) , (12) 

we get in our metrical spacetime the following results: 



X) = ,l^Yi^m{e,v)e-"^^J,Xpr) (13) 
and 

= aV + , (14) 
where Yim{0,(f)) are the spherical harmonics, J^^ the Bessel functions and 



ui = a~^^{l + 1/2)2 ^ 2(1 - a2)(^ - 1/8). So, according to (|0) our zero 
temperature heat kernel is given by 



I TT [ c/p^0<x(a;)0<x(a;')e 
J-oo 2n Jo TZ 



l,m 

2 I ,'2 



g-(Ar)V4s e-^i^ - ( rr' \ 

(15) 



i=0 



where we have used the addition theorem for the spherical harmonics, with 7 
satisfying the well known relation with the original angles (^, </)) and (^', (/)'). 
In (^) ly^ is the modified Bessel function. 
Comparing (|15D with (||) we can identify 
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and 



Now we want to call attention for the fact that for a = 1, ui become equal 
to / + 1/2 and it is possible to get a closed expression for the sum in (0) 
[pj]l , getting the ordinary heat kernel function for the flat-space case 



1 p-{x-x'f/4s 

^L"^^H-'-'^) = I^^^^- (18) 



Now we are in position to calculate the heat kernel function at non-zero 
temperature in the global monopole spacetime. According to only 
ii'(i)(r, r'; s) will be affected by temperature, so we get 

-| oo 

KiMr, r'; s) = -j= e-(---'^+"/^)V- . (19) 

Combining Eqs. ([1^) and (p!?!), we have Kf^^x, x'\ s) and consequently our 
thermal Green function, (^/^(x,^'), using (^). Our result is: 



-j^ oo 

Gp{x,x') = GUx,x') + J2 E(2^ + l)Qu,-i/2iufs)Piicos^) , (20) 

avr rr „_^oz=o 

where 

a2(r-r' + ra/3)2 + r2-hr'2 , . 

«/3 = > (21) 

and 

-| oo 

GUx,x') = -^^(2/ + 1)Q,,_i/2(moo)PKcos7) ■ (22) 

(Unfortunately it is not possible to write this thermal Green function in terms 
of a single special function.) 

Our temperature independent Green function Goo{x,x') coincides, up to 
a redefinition of the "Euclidean" time r by ar, with the Green function for 
a massless scalar field at zero temperature given in Ref. B. 
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For the case a = 1, Eq. (|20| ) provides 

^ ^ ' 47r/?|f -F| cosh (^|f-F|) -cos (^(r-rO) 

From the expression above it is possible to obtain its zero-temperature hmit, 
taking (3 — * oo. In this hmit we get the ordinary Euchdean Green function 
given in the hterature: 

Gix, x') = A77 — ^—T^ . (24) 

^ ' ^ 47r2 (x - x')2 ^ ^ 

3. The Computation of {(f)^(x))i3 at Nonzero 
Temperature 



The thermal average of {(j)'^{x))i3 can be obtained computing the coinci- 
dence limit of the Euclidean thermal Green function as shown below, 

(0^(x))/3 = lim Gp{x,x') = lim [Goo(a^5 a^') + Gf^^x, x')], (25) 

where we have separated the purely thermal part of (pOD , defined as Gjsi^x, x'). 

The above procedure gives us a divergent result which comes exclu- 
sively from the zero-temperature contribuition of the thermal Green func- 
tion, Goo{XiX), which is proportional to the behavior of the Legendre func- 
tion evaluated at the unity. So, in order to obtain a well defined value for 
{(j)^{x))f3, we extended to our thermal problem the renormalization procedure 



given in Ref. [0: we subtract from the complete thermal Green function 
the Hadamard function, G^^\x,x'), which is given in terms of the square of 
the biscalar geodesic interval, a{x, x'). So, the renormalized thermal average 
value for {(l?{x))p is given by 

{<p\x))p^R,r.. = \vm[Gp{x,x')-G^^\x,x')] , (26) 

x'^x 

which is explicitly equal to 

1 00 00 / 2o2 

n=l 1=0 \ 
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In their paper, Mazzitelh and Lousto^ have computed only the zero 
temperature contribution, {4>'^{x))oo,Ren., up to first order in the parameter 

= 1 — a^, considered smaller than unity. So, we will not repeat their 
calculation. In this paper we compute the thermal contribution for the square 
of the scalar field operator up to the first order in the parameter rj'^ like inP]: 

1 oo oo / 772/92 \ 

= ^ E E(2/ + l)Q.-i/2 [I + j • (28) 

In order to do that let us take the integral representation for the Legendre 



function 13 



1 c~'^^ 

Q„_.M00Shrt = ^/^ *^=g=g=. (29) 

Substituting this expression in (|28|) , it is possible to proceed the summation 
on the angular quantum number, by using the expansion for i>i, 

After some other intermediate steps we obtain 



where 



and 



= ^(1 + v') - + ,^^s,, (30) 



Si = 12 dt , , = , ' , (31) 

r^^Jpri vcosht — coshor, smh(t/2) 



00 /-oo 1 + 

S2 = Y. dt ^ ^ I (32) 

^i^pn V cosh t — cosh p„ sinh (t/2) 

with the notation p„ = arccosh(l + n^/9^/2r^). 

Unfortunately it was not possible to obtain an explicit result for both 
integrals above in terms of elementary or special functions. The best that we 
could do was to proceed a numerical evaluation of both integrals for specific 



^In the Appendix A we present in more details the steps adopted by us in this 
calculation. 
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values of the ratio ( := (3/r and different values of n. In our numerical anal- 
yses the results were developed in the high temperature, or large distance 
limit, C -C 1. They are shown in Figs. (la) and (lb) for Si{() and 5'2(C) 
respectively. From the graphs displayed in Fig. (2), which exhibit the loga- 
rithmic behavior for the previous figures, it is possible to infer the following 
dependences for Si{() and S2{() with (: 

SiiO = ^ (33) 

and 

^^(0 = ^ , (34) 

where ci = 14.26 with gi = 0.92 ± 0.07 and ca = 17.82 with ga = 2.02 ± 
0.02. From these results we can see that the most relevant contribution to 

2 

{(f> )/3 is given by the term independent on parameter ^. It is, approximately, 
proportional to 1//?^ and consequently independent of the distance from the 
point to the global monopole. 



4. The Thermal Average of {Tqq(x 



The energy-momentum tensor, T^^y{x), is a bilinear function of the fields, 
so we can evaluate its vacuum expectation value, (T^,y(a:))/3, by the stan- 
dard method using the Green's function^]. The thermal vacuum average. 



consequently, can also be obtained using the thermal Green function. 

In Ref. the general structure of the renormalized VEV of the energy- 
momentum tensor associated with a massless scalar field in the global monopole 
spacetime at zero temperature is presented. (There, it was obtained consid- 
ering the VEV of square of the field operator on the base of dimensional 
arguments, symmetries and trace anomaly.) The objective of this section is 
to calculate the nonzero temperature correction to this vacuum polarization 
effect. So we concentrate on the thermal correction. The purely thermal 
average of energy-momentum tensor is given by: 
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lim 



(1 - 20 V^V..G^(x, x') - 2^V,V.G/,{x, x') 



+ [2^ - 2j 9f.u{x)9''' {x,x')VpV^'Gf3{x,x') 
-iG^y{x)Gp{x,x') - 2^^g^i,R{x)Gp{x,x'] 



(35) 



where G^^^ and R are the Einstein tensor and the scalar curvature, respec- 
tively. 

As we have already mentioned, we shall consider, for the sake of simplicity 
only, the zero-zero component of (|35D; so, taking fi = u = and using the 
fact that dt'G{x,x') = —dtG{x,x'), we obtain 



lim 

x'^x 



-dlGp + (2e - ^) gA^)g'"'\x,x')VpV.,Gp{x,x') 



+-ai-mR{x)Gp{x,x') 



(36) 



On the other hand we know from ([2^) that 
Gp{x,x') 



7 E E(2^+l)^?^.-i/2 ^ 7^, i^Kcos7). 

(37) 



So, after some intermediate steps, we arrive at 

2 oo oo 



(roo(a;))/3 



^^(2/ + l)[S„(a,C)Q 



n=l i=0 



(2) , 
ui-l/2\ 



n(«, e, C)Ql;^/2(^n) + ^n(«, (^n)] , (38) 



where 



Sn(a,C) 
'^'n(a,C,C) 



r 2 2.2^Z\[2C(A'-4)--(r^V + 4)] 



1 



nCa^^{n^C + 4) 
^ '(Sn^C^ - 2)a^ 



{2e [(Sn^C^ - 2)0" - 2 



(40) 
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and 



Here 



*„(a,{)= 2? 



(41) 



1 + and C = In Eq. ([ 



(^^'"^(-z) are the associated 



Legendre functions given by 



for integer positive m. 

The equations (^H|)-(^) give a formal expression for the purely thermal 
correction to vacuum expectation value of (Too(x));3. However from these ex- 
pressions it is not possible to obtain a concrete conclusion about its behavior 
with temperature. So, in order to provide this quantitative information, let 
us expand (Tqo)/? in powers of the parameter rf up to the first order. How- 
ever, the dependence of (Too(a;))/3 on = 1 — appears in two different 
ways: (i) in the coefficients (^)-(|4l|) and (ii) in the argument of the Legen- 
dre functions. After long calculations^ which also include the summation in 
the angular quantum number /, we obtain an expression containing a large 
number of simple algebraic functions and integrals similar to the previous 



ones given in (31) and (32). Our final result is 



(Too(x))^ 



'30/3^ ^ V24 ^Qj r2/32 



+ 



TT 



11 1\ 1_ 

180 ^^9; /34 + 



(42) 



where 



^(0 = ^(1-27^-^^ 



(43) 



The function -F(^, () is expressed in terms of integrals as we show below 



2^,4 



4 



/o,3(c) + e/o,i(c) - e/2,3(c) - eV2,i(c) 



^3(C) + e^i(C) 



2^ 



^/^4,1,5(C)-4^1A3(C) 



(44) 



^In the Appendix B we present these calculations in more details. 
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where 



/m(C) = E(3^'c' 



n=l 



d 



udu 



1 



n=l 



9 AO 

and 



Ci2 



rfy yarccosh(y) Vcosh u-y smh!' {u/2) J j^=i+„2(^2 /2 
^ urfw 1 



d?/^ 7arccosh(y) i/cosli m — ?/ sinh-' {u/2) 



(45) 
(46) 



E 

n=l 



/2 + pe"" 



arccosh(y) a/ cosh u — ?/ sinh'^(u/2) 



(47) 



1, 3; / = 1, 2; p = 0, 1 and q 

t2a2 a\ 1/ 



3,5. We call attention for 

■2 



with i = 0,2; j 

the fact that = 2^{nX^ - 4) - ^n^C^ + 4) and = ^2(2/2. 

Unfortunately, we could not express the integrals above in terms of ele- 
mentary functions. Once more, in order to provide a quantitative information 
about the behavior of (^2]), we have to proceed a numerical analysis of the 
integrals as a function of ( = (3/r. However, the numerical calculations be- 
come more complicated here because of the derivatives which come from the 
development of Ql^^ and Q^J'^ ■ So making an appropriate transformation of 
variable, u := arccosh (^^^ + 1^, the new integrals in the variable q present 
fixed limits. This procedure allows us to take their derivative in a simple 
way. In the high temperature regime, C <^ 1, we have obtained the following 
numerical results: 



30/?4 ^ (24 



6 / r2/32 



+ 



^2(0 



where 



TT 



^i(0 = -Y^(ll + 520 , 



^^2(0 



768 



(64^2 ^ 1^25^ - 832) 



(48) 

(49) 
(50) 



and 



4V2' 



2^,4 



TT^r 



^/o,3(c) + e/o,i(C) - e/2,3(c) - e72,i(C) 



12 



+ 2^ 



^/i2,l,5(C) - 1^1,0,3(0 



(51) 



Here the behaviors of all the above expressions with ( = jS/r are given below 

n 

(i) For the functions QjiC) we have: 
~9i = -^ 



132.71 10.19 8.11 
+ — + 



28.20 2.26 8.11 
+ -:— + + (52) 



and 



^3 



,129.84 10.19 
<—:— + 



(53) 



with ai = 2.26 ±0.04, bi = 2.03 ±0.03, ci = 3.00 ±0.00 and 02 = 2.22 ±0.02, 
62 = 1.96 ± 0.23, C2 = 2.99 ± 0.01. Also 03 = 4.01 ± 0.01 and 63 = 4.0 ± 0.00. 
(ii) For the functions fij we have: 



/( 



55.21 



03 



a4 



with 04 = 1.96 ± 0.01, 



21 



118.85 21.27 

± 



with as = 2.27 ± 0.07 and 65 = 3.04 ± 0.03, 



48.87 915.95 



23 



with ae = 3.99 ± 0.04 and fog = 2.13 ± 0.07. 
{in) Finally, for the functions hi^p^q{(), we have: 

77.63 



(54) 



(55) 



(56) 



^2,1 



and 



,1,5 



0,3 



16.36 



(57) 



(58) 



^We exhibit only the contributions which depend on C with power higher than or the 
same order as 0{1/C^). 
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where ay = 2.00 ± 0.02 and 67 = 3.97 ± 0.03. 

From our numerical results, which were found for C, = jS/r belonging 
to the interval [0.01,0.1], we can conclude that most relevant terms are of 
order 1/C^- According to (|51|) these terms represent contributions to Eq. 
( ^2]) which are proportional to 1/13'^ and are of the same order of the leading 
term, obtained analytically and independent of the distance. So the thermal 
average {Tqq{x))p can be written in the form below, 

{TM)l={T,,{x))l + T'g{^\ (59) 

where g is a. function of the dimensionless variable 7 = rT, whose expression 
has been found by us approximately up to first order in the parameter rf. 

5. Concluding Remarks 

In the present paper we have obtained the Euclidean thermal Green func- 
tion Gp{x^x') associated with a massless scalar field in a global monopole 
spacetime. Our result was expressed as a sum of a zero-temperature Green 
function plus a purely thermal contribution. Using this Green function we 
were able to obtain the corrections due to temperature in (0^(x))/3 and 
(Too(a;))/3. Because in both calculations, the thermal contributions could 
not be expressed in a simple way, but some of them in terms of not solvable 
integrals, their dependence on temperature could not be evaluated explicitly. 
So, in order to get some quantitative information, we had to develop a numer- 
ical analysis. We decided to consider these analyses for a high-temperature 
regime. Our final results for the thermal vacuum average {(p'^^x))^ and (Tqo)/?, 
present the following behavior: 

mx))l={<p\x))l + T'f{^) (60) 

and 

{no{x))l= {no{x))l + T'g{^). (61) 

For both case the functions / and g, which depend on the specific geom- 
etry, were found by us up to the first order in the parameter rf = 1 — a^, 
considered smaller than unity. In fact for a typical grand unified theory the 
parameter rjo is of order lO^^GeV, so r/^ = Sttz/q ^ 10~^. 
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A few years ago, Smith|T6| and Linet|T^ have obtained the Euchdean 
thermal Green function and also the thermal average (0^(x))/3 and (Too(x))/3, 
for a massless scalar field in a cosmic string spacetime. For this case, Linet 
could present his results in a closed form. These two beautiful papers give us 
the motivation to analyse similar phenomena in a global monopole spacetime. 

We want to emphasize that in all the above numerical evaluations, we 
considered sufficiently large number of terms in the summation in order to 
obtain good approximate results. (Independent numerical simulations pro- 
vided us the minimal values for n in each series.) 

The obtained results may have some applications in the early cosmology, 
where the temperature of the Universe was really high. We can see that 
the thermal contributions to and (Too(x))/3 modify significantly the 

zero temperature quantitaties in that epoch. Because of this, they should be 
taken into account, e.g., in the theory of structure formation. 
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A High-temperature limit of the thermal av- 
erage of 



In this appendix we give a brief explanation of our procedure used to 
obtain the high-temperature limit of the thermal average of Ac- 
cording to (p7|), the purely thermal contribution to this quantity is: 



n=l 1=0 



with 



1 + 



2/-2 



and 



A 



(A.62) 



(A.63) 



2 2 ' 

where ( = f3/r. 

Expanding first the argument of the Legendre function in ([A.62| ) up to 
the first order in 77^, we get: 



n=l 1=0 



d 



{AM) 



We can substitute z// = 1 + 1/2 into the second term of (|A.64|) . This 
allows us to make the summation in the angular quantum number I and 
consequently make the second summation in n, resulting 



4^2j.2 



(A.65) 



12/3^ ' n=U=0 

Using the integral representation (pQ]) for the Legendre function and the 
expansion 

1\ /- f] 



1 + 



1 



2e - 1/4 2^ 
2/ + 1 ' 
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we can obtain a simpler expression for ( |A.65| ) by the summation in / as shown 
below: 



yf Ing 

(W) 



(A.66) 



where q = e~*. (This summation was also present in Ref.|^. There is, 
however a misprint, the factor rj'^ appears dividing by 2, for this reason we 
reproduce the correct result. We want to call attention that this misprint, 
has been corrected in further calculations and no mistake was found by us). 
Now going back to (|A.64|) we get the following expression 



2 

(0 



[i + v' 



:S2{0, (A.67) 



12/52' ' ' ' 87r2r2 "'^^ ' 327r2r2 
where Si{() and 5'2(C) are integrals given by (|3lD and (|3^) respectively. 
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The expansion of (|SS| ) in powers of 77" 



The expansion of (|38D in powers of r]'^ is very long. In order to make our 
explanation clear to the reader, let us sepatate this calculations into three 
parts: 

(2,) 

(i) The first one refers to the term in Ql/_i/2i 



;i - E E(2/ + l)Sn(l - n\ e, C)Q2^v2(^n - ri'Xl). (B.68) 

n=l 1=0 



Making the expansion in the argument of the Legendre function in powers 
of rf, we get: 



n=l 1=0 



n^Q + 4 



dy 







y=z„_ 





where 



S„(x,C)=2e(n^C'-4)--(nV + 4) . 



[B.69) 



(B.70) 
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(2) 

Using the definition for Q] and the integral representation for the Leg- 
endre function, we obtain, after the expansion in powers of rf and the sum- 
mation on /, the following expression: 



h 



vl(i + 40 



-TT 



30 "^^45/ 13^ 



rf + 



(B.71) 



where the functions 5'i(C) and 5'2(C)(-2n) are given in (^61). 

The other two terms of (|38D are developed in similar way. After some 
intermediate steps we obtain the following results: 
{ii) For the term proportional to Q^l^_xi2 have found 
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(B.72) 



where the functions /o,3(C), /o,i(C), /2,3(C) and /i,3(C) are given in Eq. 
(m) For the term proportional to Qui-i/2{zn) we have: 
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(B.73) 



where the functions /i2,i,5(C) and /ii,o,3(C) are given in the Eq. 
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Figure 2: This figure exibit the logaritmic behavior for both functions, 5*1 
and 5*2 with (, in the region [0.01, 0.1]. From it, it is possible to estimate the 
leading dependence for both quantities with 
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